Magneto-electric coupling in zigzag graphene nanoribbons 
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Zigzag graphene nanoribbons can have magnetic ground states with ferromagnetic, antiferromagnetic, or 
canted configurations, depending on carrier density. We show that an electric field directed across the ribbon 
alters the magnetic state, favoring antiferromagnetic configurations. This property can be used to prepare ribbons 
with a prescribed spin-orientation on a given edge. 



Introduction — Expanding techniques that can achieve elec- 
trical control of spin is a key goal of both metal and semi- 
conductor spintronics.Hl 0] In metal spintronics, for exam- 
ple, electrical spin-transfer torquel]! research seeks to am- 
plify the potential|4] of technologies based on giant magne- 
toresistance Q\ and tunnel magnetoresistance. J@] The aim 
of research on dilute magnetic semiconductors|7|] is to cre- 
ate semiconductor materials in which magnetic properties are 
as sensitive to doping and external gate potentials as electri- 
cal properties. Recent interest in the spin Hall effect |@] and 
the topological magnetoelectric effect |@1 is motivated by a 
search for effects which enable electrical control of spin in 
non-magnetic materials. In this context it is interesting to ad- 
dress the possibility of interesting magnetoelectric effects in 
graphitic material. The physics of zigzag graphene ribbons 
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Magnetism is expected in any graphitic material containing 
ribbon segments with zi gzag edge | flfl [III [H [H H 
13, GJ, El HI EH |22L M uA |25literminations, for exam- 
ple highly defected bulk graphitic material. 14 ll l42l |43L 14411 
Thanks to progress in structural control of graphene flakes and 
related materials [45], prospects for mastering graphitic mag- 
netism have improved. The perfect zigzag nano-ribbon stud- 
ied in this Letter may be viewed as a model system in which 
graphitic magnetism is exhibited in its simplest and most es- 
sential form. 



The magnetic ground state of a zigzag nanoribbon has col- 
lective moments localized near its edges. In the absence of 
a transverse electric field, a doped nanoribbon has either full 
or partial orientational alignment ll46l 14711 between moments 
on opposite edges, as illustrated schematically in Fig.[l]. We 
show that an external electric field applied across the ribbon 
can control the relative orientation angle 0, and that this prop- 
erty can be used to prepare ribbons with a prescribed spin- 
orientation on a given edge. 

Hubbard Model Mean-Field Theory — Since ab initio density- 
functional and simpler model Hamiltonian approaches make 
essentially identical predictions IU9U231 l24l 14311 . we base our 
analysis of zigzag ribbon magnetism on a Hubbard model 
which allows the underlying physics to be identified more 
clearly. The success of the Hubbard model has been shown 
to be due to the essentially local character of edge magnetism 
in graphene ribbons [51]. The Hubbard model mean-field 




FIG. 1: Schematic representation of the influence of a lateral elec- 
tric field on zigzag ribbon magnetism. Upper panel. Non-collinear 
ground state of a doped zigzag ribbon illustrating the angle 9 be- 
tween moment orientations on opposite edges. The red box indicates 
the unit cell of the ribbon. Lower panel. A lateral electric field drives 
the angle 6 to %, an antiferromagnetic configuration similar to that 
of an undoped ribbon. 



Hamiltonian, 
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has a term which represents hopping between nearest neigh- 
bor 7T-orbitals with amplitude 70 = 2.6eV, an external poten- 
tial term which accounts for the transverse electric field, and 
a mean-field interaction term. The operator cj a creates a n or- 
bital electron at site i with spin a, S is the transversal electric 
field, yt is the position of lattice site i along the ribbon width 
and x a i a represents the elements of the three Pauli matrixes. 
As a convenience we include a constant term v ext = — U in 
the external potential which removes the interaction with a 
unit charge on each site from the mean-field quasiparticle en- 
ergy. (All spin indices in Eq. (T) are summed over.) Note 
that the mean-field interaction energy of an electron on site 
i is spin-dependent and proportional to the density of oppo- 
site spin electrons. Following Yazyev et ai, [43] we choose 
U = 3eV a value slightly larger than estimates based on the 
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local density approximation used in some previous Hubbard 
model analyses IJjl 21, 47]. We have used 1200 fc-points in 
the Brillouin zone for the self-consistent calculations. Zig-zag 
edge magnetism is very sensitive ll47ll to the net charge density 
of the ribbon 8n, which we refer to as doping whether due to 
chemical dopants or gate voltages and measure per repeat dis- 
tance a = 2.464 along the edge. The corresponding areal den- 
sity is 8ri2D = Sn/W where the ribbon width W = V3Na/2 
and N is the number of atom pairs per ribbon unit cell. 
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FIG. 2: Edge state bands of a zigzag ribbon with N = 20, nearest- 
neighbor hopping 70 = 2.6eV and on site repulsion U = 3eV under 
a transverse electric field of S = 0.2eV /nm. Left panel: In the AF 
case high-spin (see text) bands shift toward the Fermi energy and 
the low-spin bands move away. The wavefunctions of the occupied 
bands shift density to the low-energy edge. Right panel: In the F case 
Ag shifts the wavevector at which majority and minority spin bands 
cross to slightly larger values of |Jk|. The dotted lines show the bands 
at Ag = 0. 



Edge-Only Model — The influence of a transverse electric field 
on electronic structure in neutral ribbons has been studied pre- 
viously for both armchair [48] and zigzag cases [ 16, 49, 50] 
using density functional theory. The essentials of zigzag- 
ribbon magnetism and of the transverse-field magnetoelec- 
tric effect are captured by an edge-state-only model lll9l 12111 : 
the qualitative discussion below refers mainly to this model 
and to the special case of collinear magnetic states (9 = or 
= k), but the numerical calculations and the phase-diagram 
results are based on solutions of the the full non-collinear %- 
band Hubbard model self-consistent field equations. In the 
edge-only model the spin-dependent mean-field Hamiltonian 
for collinear states takes the form ll2ll 51 ] 



H a (k) = (aA z (k) + Ag (k))T z + (A (jfc) + h z )al + 1 (jfc) i x . 

(2) 

Here the x a are Pauli matrices which act on the which edge de- 
gree of freedom. The terms proportional to x z in Eq. J3 there- 
fore represent the difference in energy between left (t z — > 1) 
and right — > — 1) edges for electrons of spin a, whereas the 
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FIG. 3: Inter-edge relative spin orientation angle 6 evaluated on the 
outermost edge atoms as a function of transverse electric field Ag for 
a series of doping Sn values. For Ag = doping leads immediately 
to a canting angle 9 < n and eventually for S n larger than ~ 0.05 to 
a ferromagnetic state with 9 = 0. A finite Ag favors the AF state as 
explained in the text and drives 9 toward its undoped value. 



term proportional to % x represents the momentum-dependent 
inter-edge hopping amplitude. Zig-zag edge magnetism fol- 
lows from the property M2 ill that t (k) vanishes rapidly with rib- 
bon width in the part of the Brillouin-zone {2n/3a < \k\ < 
n/a) in which edge states reside. In Eq. (O, A z (k) captures 
the difference between exchange energies on opposite edges, 
which vanishes in the 9 = (F) state in the absence of a trans- 
verse field, whereas Ao (k) captures the spin-dependence of 
the edge average, which vanishes in the — % (AF) state. 
(An irrelevant spin and edge independent exchange energy 
has been dropped from H a .) Both exchange energies are large 
only for the edge-states (\k\ > 2n/3a ). h z accounts for Zee- 
man coupling to the ribbon spins by an external magnetic field 
when present. 

Magnetoelectric Coupling in Undoped Ribbons — The 
eigenenergies of this Hamiltonian are 



E± (jfc) = ct(A (Jt) + h z ) ± yj (cjA, (Jfc) + Ag {k)Y + f 2 (Jt). (3) 

Note that there are always four distinct eigenvalues in the fer- 
romagnetic case, whereas the antiferromagnetic state bands 
occur in doubly-degenerate pairs when Ag — > 0. For undoped 
ribbons the lowest two edge states bands are normally fully 
occupied. In Fig. 2 we plot ribbon band structures for both 
AF and F states of a neutral ribbon calculated using a con- 
stant transverse electric field of $ — Q.2V/ nm. States that are 
shifted down (up) in energy relative to the Ag = case are 
localized on the low (high) potential side of the ribbon. (Note 
that a constant field generates a ^-dependent Ag because of 
the ^-dependence of the degree of edge state localization.) In 
the F-state the Fermi energy is pinned to a band-crossing near 
\k\ = 2n/3a between the higher energy majority spin band and 
the lower energy minority spin band. In the AF case we refer 
to the spin-orientation which dominates occupied states on the 
low (high) potential side of the ribbon as the low-spin (high- 
spin). (In Fig. 2, { is the low-spin.) A transverse spin shifts 
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FIG. 4: Energy differences (per ribbon unit cell) between the spin 
collinear AF, F solutions and, for finite-doping, the minimum energy 
non-collinear (NC) solution. The energy difference between AF and 
F solutions for Sn = has a weak transverse electric field depen- 
dence. 



the energies of both occupied and unoccupied high-spins to- 
ward the Fermi level, lowering the gap. A sufficiently large 
transverse field will close the indirect gap, creating a half- 
metallic band structure with only high-spin bands crossing 
the Fermi level. This is the magnetoelectric effect discussed 
m earlierOil IIIE}] work. The energy difference between 
F and AF states 12 ill is relatively unchanged by a transverse 
field. Below we show that in doped ribbons a transverse field 
tilts the competition between F and AF states in favor of the 
latter, yielding a distinct and stronger magnetoelectric effect. 
Magneto-electric coupling at finite doping — We consider for 
definiteness the case of n-type ribbons in which carriers are 
added by gate doping. In the absence of a transverse electric 
field, doping favors| 4& 47 ] the gapless F state over the gapful 
AF state. The transition between undoped AF and doped F 
states occurs continuously by varying the relative orientation 
angle 9 between its AF (9 = it) and F (9 = 0) end points. The 
effect we discuss in this paper is based on the following simple 
observation concerning the edge-state bands plotted in Fig. 2. 
In the AF case the conduction band states which are occupied 
upon doping are high-spin antibonding states, which are local- 
ized on the low-energy side of the ribbon. For the F state, on 
the other hand, there are occupied states in two bands, one lo- 
calized on the high-energy side and one localized on the low- 
energy side. The net effect is that a transverse field favors the 
AF state in doped zigzag ribbons. In Fig. 3 we plot the relative 
angle between spin polarizations on opposite edges vs. lateral 
electric field for a series of different doping values. These re- 
sults were obtained by non-collinear spin self-consistent field 
calculations and confirm the expected magnetoelectric effect. 
In Fig. 4 we compare the transverse field dependence of the 
energy difference between F and AF states for doped and un- 
doped systems. The electric field strength required to convert 
F states into AF states in doped ribbons is much smaller than 
the field required to close the AF-state gap in the undoped 
case. 

Discussion — Typical results for the edge-state bands of both 
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FIG. 5: Top panel: Band structure of n and p gate-doped N = 8 
zigzag graphene nanoribbons in the presence of a transverse electric 
field, $ = 0.2 V/nm, strong enough to produce a collinear 6 — % 
ground state. Note that the edge-state bands are half-metallic in both 
cases. Bottom panel: Spin resolved occupation n (y) across the 
robbon: majority spin electrons (<J =f for n-doped and cr =| for p- 
doped) in this figure accumulate on the high potential edge of the 
ribbon for Sn = 0.06 and on the low-potential side for Sn = —0.06. 



electron and hole doped zigzag ribbons with a transverse field 
strong enough to induce the 9 — % state are illustrated in 
Fig. 5. Because of the partial occupation of the highest un- 
occupied band of the n-doped case and the lowest unoccupied 
band in the /?-doped case, the 9 = n state has an overall spin- 
polarization proportional in magnitude to the doping. Note 
that these 9 = n states are always half-metallic. This fer- 
romagnetic component of the order allows Zeeman coupling 
from an external magnetic field to fix the spin-orientation on 
each edge. It follows from Fig. 5 that the majority spins are 
high-energy spins in the «-doped case and low-energy spins 
in the p-doped case. For a fixed magnetic field direction, the 
spin-orientations on both edges can therefore be switched with 
a gate voltage, which changes the sign of the carrier density. 
This remarkable property of zigzag edge magnetism has no 
parallel of which we are aware in any other magnetic system. 

In closing we remark that the one-dimensional character 
of zigzag edge magnetism works against robust collective 
spin properties. As analyzed in more detail elsewhere lfl8i 
I22I1 . the consequences of reduced dimensionality are some- 
what mitigated by the the substantial stiffness lfTii H2l of 
zigzag edge moments. Nevertheless, robust magnetism in 
graphitic nanostructures will likely require exchange-coupled 
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two-dimensional ribbon networks. The rather unique proper- 
ties of graphitic magnetism discussed in this Letter motivate 
an effort to realize structures of this type. 
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